Particle Entanglement in Rotating Gases 
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In this paper, we investigate the particle entanglement in 2D weakly-interacting rotating Bose 
and Fermi gases. We hnd that both particle localization and vortex localization can be indicated 
by particle entanglement. We also use particle entanglement to show the occurrence of edge re- 
construction of rotating fermions. The different properties of condensate phase and vortex liquid 
phase of bosons can be reflected by particle entanglement and in vortex liquid phase we construct 
the same trial wave function with that in [Phys. Rev. Lett. 87, 120405 (2001)] from the viewpoint 
of entanglement to relate the ground state with quantum Hall state. Finally, the relation between 
particle entanglement and interaction strength is studied. 

PACS numbers: 03.67.Mn, 03.75.Gg, 73.43.-f 



I. INTRODUCTION 

Quantum entanglement, which is considered to be the 
most non-classical phenomenon in the quantum world, 
has been attracting much attention in last decades. On 
the one hand, it has been identified as a key resource 
in many aspects of quantum information theory, such 
as quantum teleportation, quantum key distribution and 
quantum computation On the other hand, the con- 
cept of entanglement is proven to be useful in condensed 
matter systems Q, such as in spin chain, Bose-Einstein 
condensate and electronic quantum Hall effect. 

As a typical condensed matter system, the two- 
dimensional (2D) rotating atom gas is a very interesting 
field because it opens a door to many important quantum 
phenomena 0, H| ■ Many excellent experimental and the- 
oreticalpapers focus on the formation and melt of vortex 
lattice [5| , many-body energy spectrum [y, 0] , its analogy 
to quantum Hall effect of electrons in a magnetic field [a- 
[i~3j . particle localization and vortex localization fl4T4l7) ]. 
the comparison between the results of exact quantum 
numerical solution and mean-field theory [l8| and so on. 
However, the entanglement in this system has not been 
investigated as extensively as in other condensed matter 
systems, for example in spin chain systems. 

In one of our previous work, we use some quan- 
tum information concepts, such as fidelity susceptibil- 
ity and particle entanglement to study the properties 
of 2D weakly-interacting rotating Bose-Einstein conden- 
sate trapped in a harmonic potential We find that 
the single-particle entanglement, defined as the von Neu- 
mann entropy of the single-particle reduced density op- 
erator, can indicate some angular momentum of real 
ground states and reflect the properties of condensed 
phase (single- vortex state) and vortex liquid phase. This 
indicates that in 2D rotating atom (both boson and 
fermion) gas, some other important properties of this ro- 
tating system may be shown by entanglement. 



In this paper, we mainly investigate the particle en- 
tanglement of 2D weakly-interacting rotating atom gas 
trapped in a harmonic potential in disk geometry. In or- 
der to compare the results of bosons and fermions, we 
choose repulsive Coulomb interaction. For a fixed par- 
ticle number N, we calculate the single- particle entan- 
glement Si of the ground states in a series of subspaces 
of fixed angular momentum L. We find that both par- 
ticle localization and vortex localization can be reflected 
by the oscillation of Si with the increase of L. This 
oscillation is universal for bosons and fermions. Then 
we study the single-particle entanglement of the ground 
states in subspaces whose angular momentum are special 
functions of N. For the subspace L = N, through the 
comparison between the results of bosons and fermions, 
we find that the edge reconstruction of rotating fermions 
can be indicated by the single-particle entanglement. For 
the subspace L _= (N - N)(N + N ~ k)/k with k > 1 
an integer and N the smallest non-negative integer mak- 
ing N — N to be divisible by k, we relate its subspace 
ground state with quantum Hall state from the view- 
point of entanglement and construct the same trial wave 
function as that in Ref.Q. Finally, we study rotating 
bosons trapped in a quadratic plus quartic trap inter- 
acting with contact potential, both attractive and repul- 
sive. We find both single-particle and two-particle en- 
tanglement of the ground state in subspace L = N and 
L = N(N — 1) increase with the interaction strength for 
fixed N. However, when the interaction is attractive, 
the single-particle entanglement of the ground state in 
subspace L — N does not decay with N, verifying that 
attractive bosons do not condensate 1201. 



II. MODEL 
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We focus our attention on the 2D rotating spinless Bose 
and polarized Fermi atom gases. In rotating reference, 



2 



the Hamiltonian of N atoms is as follows: 

N h 2 N 

n = E { - + v (^) - L ^ Q ) + E w ( r i» r i)= 



i<j=l 



(1) 



where V(r) is the trap potential, — L z f2 is the rotating 
energy with rotation frequency Q and £Y is the interac- 
tion energy. U > (<)0 represents repulsive (attractive) 



interaction. Typically, V(r) 



,2„2 



is a harmonic po- 



tential. U can be either a short-range contact interaction 
S (r j — r j ) or a long-range Coulomb interaction 1/ |r j — Tj | . 

In the weakly-interacting limit, we can use single- 
particle lowest Landau level (LLL) wave function (fi(z) — 
_L^ z ' e -|z| /2 ; wnere z = (x + iy)yJrmo/h is the di- 
mensionless position of a particle, to make second- 
quantization of Eq.([T]). The obtained Hamiltonian in the 
subspace of fixed angular momentum L z = Lh takes the 
form: 



n L = (L + N)hiu - Lhtt + U 



L 

E 

i,j,k,l=0 



Ui J j,k,ia,i a i a k(''i,(2) 



where \Uq\ represents the strength of interaction and 



Uij,k,i = jj dz 1 dz2ip*(z 1 )ip*(z2)U(z ll z 2 )ipk(z2)v>i{zi). 

For contact interaction, the expression is simple: 

Ui,j,k,l = wfj S^wF. §i +i' k + l - But for Coulomb mter " 
action, the expression is a little complicated [2l|] : 



U, 



= 8. 



i+j,k+l 




ilkl T(i + j + 3/2) 



where 



A, 



B, 



r\ r(i + l/2)r(i + t + l/2) 

^ »!(r - i)\ (i + t)W(i + s + 1 + 3/2) ! 

H r(z + i/2)r(i + t + i/2) 

i\(r - i)\ (i + ty.T(i + s + 1 + 3/2) 



E 



E 



x(2i + * + l/2). 



r(x) = J +o ° t x_1 exp(— t)dt is the usual Gamma func- 
tion. 

After diagonalizing Eq.Q numerically by Lanczos al- 
gorithm, we obtain its subspace ground state in 
angular momentum Fock representation. Here we con- 
sider entanglement between particles [22j (In Ref.[23j], 
the entanglement between angular momentum orbits is 
investigated to obtain the topological entropy in rotating 
Bose-Einstein condensate). We define the single-particle 
reduced density operator of \§>l) as (pi)ij — ^(a^aijL — 
jj§i,j{a\ai) L , with {-) L = (<S> L \ ■ \$ L ). The single-particle 



entanglement of is just the von- Neumann entropy 
of p 1 : 
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Si = -Tr(pi lnpi) =lniV- y ]{a\ai) L \n{a\ai) L . 

For bosons, it's possible that (ajai)^ = N for one i and 
(ataj)x, = for j ^ i. Therefore Si > 0. For fermions, 

because {a\ai)L < 1, we can find that Si > In AT. More- 
over, considering the dimension of single-particle Hilbert 
space is L, we have 5*1 < InL for both bosons and 
fermions. 

Similarly, we can define the two-particle re- 
duced density operator of a s {p2)ij,kl = 
n(n-i) ( a k a l a 3 a i) L$i+j,k+l an d its two-particle en- 
tanglement as S*2 = — Tr(p2lnp2)- 



III. SINGLE-PARTICLE ENTANGLEMENT 

In this section we consider single-particle entanglement 
of both bosons and fermions interacting with each other 
by repulsive Coulomb potential. 

First we fix the particle number N to calculate the 
single-particle entanglement in a series of subspace 
ground states It's known that at extreme angular 

momentum L, a few particles {N is small) in a harmonic 
trap localize to Wigner molecules, independent of their 
bosonic or fermionic statistics. While at moderate 
angular momentum, if the particle number N is much 
larger than the number of vortices, the localization 
happens in vortices. The localization of particles (vor- 
tices) can be studied by particle (hole) pair correlation 
function and can be reflected by the regular oscillation 
of quantum many-body energy spectrum with L [l5|, [lj| . 
As we show below, the single-particle entanglement of 
$l) also oscillates with L and can reflect particle and 
vortex localization very well. In Figfl] we fix particle 
number N to calculate single-particle entanglement 
Si of the ground state of bosons in every angular 
momentum subspace L > 2. With the increase of L, 
Si not only has the tendency to become larger but 
also shows oscillation. To see the oscillation clearer, 
we also plot InL — Si, the difference between Si and 
its upper bound, in FigfTJ For N — 6, a series of 
local maxima of InL — Si appear: £=6,10,12,15, 
18,20,25,30,36,40,42,45,48,50,55,60,65,70,75 etc. 
The angular momentum of the real ground state is 
L g =6, 10,12, 15,20,24,30,36,40,45,50,55,60,65,70,75 etc. 
Therefore similar to the case of bosons with contact 
interaction, the local maxima of In L — Si of the subspace 
ground states can indicate some angular momentum 
of real ground state [llj]. One can note C can always 
be written as 5k or 6k with k a positive integer. This 
reflects the spatial symmetry of the particles in the 
subspace ground state, namely particle localization. The 
series C = 5k is associated with a (1,5) pentagon ring 
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structure with one particle at the center, while the series 
jC = 6k is associated with a (0,6) hexagon ring structure. 
When L is large enough (L > 50), the series C = 5k 
dominate. When L is relatively small, the series C = 5k 
and C — 6k compete with each other and the spatial 
symmetry of particles can be observed by full iV-point 
correlation function of the subspace ground state [l7| . 
For relative large N such as N — 12 and N = 20, the 
oscillations of Si and In L — Si are different from those 
for N — 6. For N — 12, Si and In L - Si oscillate 
with a period of Pl = 2 between L = 20 and L = 24 
corresponding to the two-vortex state and they oscillate 
with a period of Pl = 3 between L = 24 and L = 36 
corresponding to the three-vortex state. For N = 20, Si 
and In L — Si oscillate with a period of Pl = 2 between 
L = 32 and L — 42 corresponding to the two-vortex 
state and they oscillate with a period of Pl = 3 between 
L = 42 and L — 54 corresponding to the three-vortex 
state. This phenomenon reflects vortex localization. 

Due to Pauli exclusion principle, the smallest angu- 
lar momentum of N fcrmions is M = N(N — l)/2, 
whose subspace ground state single-particle entangle- 
ment is IniV. Therefore when considering fermions, 
sometimes it is convenient to use AL = L — M and 
A Si = Si — In AT. In Fig[5J we fix particle number N to 
calculate single-particle entanglement ASi of the ground 
state of fermions in every angular momentum subspace 
AL > 0. The oscillation of ASi is very clear so we do not 
need to plot the In L — Si as what we did for bosons. On 
the one hand, there is similarity between the oscillations 
of subspace ground state single-particle entanglement of 
bosons and fermions. For N — 6, the positions AC of 
local minima of fermionic ASi are nearly the same with 
those C for bosons, reflecting the particle localization in 
the subspace ground state. For N = 12 and N = 20, ASi 
oscillates with a period of Pl = 2, 3 and even 4 succes- 
sively, reflecting the vortex localization. On the other 
hand, there also exists difference between boson and 
fermion case. For bosons, the average angular momen- 
tum per particle ^2(3) a t which Si begins to oscillate with 
a period of Pl = 2(3) is nearly a constant. When TV = 12, 
we have £ 2 = 20/12 « 1.67 and £ 3 = 24/12 = 2.0. When 
N = 20, we have i 2 = 32/20 = 1.6 and £ 3 = 42/20 = 2.1. 
This is consistent with the conclusion of mean-field the- 
ory that two vortices form at £2 « 1.7 and three vortices 
form at £3 w 2.1 But for fermions, £2(3) is not a 
constant. When N = 12, we have i 2 = 14/12 « 1.17 
and £ 3 = 24/12 = 2.0. When N = 20, we have 
£ 2 = 20/20 = 1.0 and £ 3 = 33/20 = 1.65. Moreover, 
from Fig|31 we can see that the change of the distribu- 
tion of eigenvalues of the single-particle reduced density 
operator of bosons with the increase of L is different from 
that of fermions. For bosons, the first two (three) eigen- 
values of the single-particle reduced density operator in 
two- (three-) vortex state of bosons are always relatively 
small with the increase of L. But for fermions, the small- 
est two (three) values in the trough of eigenvalues of the 
single-particle reduced density operator in two- (three-) 




FIG. 1: (Color Online) The subspace ground state single- 
particle entanglement S\ and the difference between Si and its 
upper bound InL for 7V=6, 12,20 bosons interacting through 
repulsive Coulomb potential. The oscillations of 5*i and InL — 
Si can reflect either particle localization (for iV = 6) or vortex 
localization (for N = 12 and N = 20). 



vortex state of fermions move towards the center with 
the increase of L. 

Next we study the single-particle entanglement of 
bosons of the ground state in some special subspaces, 
whose angular momentum are functions of particle num- 
ber N. The first subspace is L = N. In Fig 01 we plot Si 
of the ground state in the subspace L = N for bosons. 
One can see that Si decays with iV monotonically, just 
like in the case of contact interaction [l9| , demonstrating 
a property of condensate phase. The second subspace is 
L = (N-N)(N + N-k)/k with k > 1 an integer, where 
N is the smallest non- negative integer making N — N to 
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FIG. 2: The subspace ground state single-particle entangle- 
ment AS 1 ! = Si — In TV for iV=6,12,20 fermions interacting 
through repulsive Coulomb potential. The oscillation of ASi 
can reflect either particle localization (for N — 6) or vortex 
localization (for N = 12 and N = 20). 



be divisible by k. For example, when k = 1, N = and 
when fc = 2, jV = for even N and TV = 1 for odd N. 
In FigEl we plot Si of the subspace ground states cor- 
responding to k = 1,2,3,4. We can see for all of them 
Si ~ ln(27V/fc — 1), demonstrating a strongly-correlated 
property of vortex liquid phase. This phenomenon im- 
plies from an entanglement view that the ground states 
in these subspaces L = (N — N)(N + N — k)/k may 
have a close relation with quantum Hall states. Re- 
calling that for a quantum Hall state in spherical ge- 
ometry, the single-particle entanglement is exactly Si — 
hx(N/v— c+1), where v is the filling factor and a is called 
shift, we can extract some information of those subspace 
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FIG. 3: (Color Online) The eigenvalues of the single-particle 
reduced density operator of bosons (left column) and fermions 
(right column). For bosons, the first two (three) eigenvalues 
of the single-particle reduced density operator in two- (three-) 
vortex state of bosons are always relatively small with the in- 
crease of L. But for fermions, the smallest two (three) values 
in the trough of eigenvalues of the single-particle reduced den- 
sity operator in two- (three-) vortex state of fermions move 
towards the center with the increase of L. For example, for 
20 fermions with AL = 26 (39), the 7th and 8th (7th, 8th 
and 9th) eigenvalues are in the bottom of the trough. 



ground states from their single-particle entanglement, al- 
though we are dealing with disk geometry. Comparing 
with our results, we can find that the ground state of 
subspace L — (N — N)(N + N — k)/k is a quantum 
Hall state with filling factor v — k/2 — limjv_ i . 00 ^^^l^ 
and when N is divisible by k the shift of its coun- 
terpart in spherical geometry is a — 2. We can ex- 
press the single-particle entanglement in another way 
Si w ln[2(iV/fc - 1) + 1]. It's known that for a Laughlin 
state of Nq particles in spherical geometry with filling 
factor v = l/m, S\ an = \n[m(N — 1) + 1]. Therefore in 
the viewpoint of entanglement, when N/k is an integer, 
our subspace ground state behaves approximately like 
a Laughlin state of N/k particles with m = 2, namely 
in our subspace ground state one particle is only en- 
tangled with N/k — 1 particles in the form of Laughlin 
state. Enlightened by this fact, we can express the sub- 
space ground state ^ k simply as a product of k Laughlin 
states, each of which consists of N/k particles and has 
a filling factor 1/2. Then considering the symmetry of 



5 



boson wave function, we obtain: (the exponential factor 
is omitted and we suppose N is divisible by k) 



N/k N/k 
i<j€Ai Km£A k 



where the set Ai(i — 1,2, fc) has N/k particles and 
the symbol S indicates symmetrization over all partitions 
of N particles into sets Ai. Now we obtain the same 
trial quantum Hall state as in Ref. Q but we achieve this 
in an entanglement manner. It has been checked that 
1 & k has large overlap with the ground state obtained by 
exact diagonalization Q . If N/k is not an integer, namely 
N ^ 0, we can slightly change the form of : 



-+i 



-+i 



i<jeAi l<meAft 

N-N N-N 

E fc 

x Yl (z P ~z q ) 2 ... Yl (z r -z s ) 2 



where the set Ai(i = 1, 2, N) has (N - N)/k + lpar- 
ticles and the set B % (i = 1, 2, k-N) has (N - N)/k 
particles. One can verify that the total angular momen- 
tum of this state is just L = N{^-)(^- + 1) + (k - 



N—N \ I N—N 



k 

1) = (N - N){N + N -k)/k. 



At last we study the single-particle entanglement of 
fermions of the ground state in some special subspaces as 
a comparison with bosons. The first subspace is AL = N. 
In Fig|4l we plot ASi of the ground state in the sub- 
space AL = N for fermions. When N is small, ASi 
decays with N monotonically similar to bosons. But 
when N > 17, AS± suddenly jumps to a higher position 
and then begins to decay again accompanied by oscil- 
lation. This is because the distribution of the eigenval- 
ues of the single-particle reduce density operator changes 
qualitatively from N = 17 to N = 18 (Fig©. When 
N < 17, the subspace ground state is a central single- 
vortex state, while when N > 17, the subspace ground 
state is a candidate of a off-center double-vortex state 
(see N = 18, 20, 22, 40 in Fig©. Therefore entanglement 
can indicate the edge reconstruction of rotating fermions 
[24| . A similar phenomenon of electrons in quantum dot 
that reconstruction of the maximum density droplet be- 
gins from the edge rather than the dot center if the elec- 
tron number exceeds N ks 15 is reported in Ref.fUJ. The 
second subspace is AL = (N-N)(N+N-k)/k. Through 
an analysis similar to that for bosons, we can obtain that 
Si f» ln[(l + 2/fc)iV — 2] and find that the ground state in 
this subspace is a quantum Hall state with filling factor 
v = k/{k+ 2) and when N / k is an integer its counterpart 
in spherical geometry has a shift a = 3 (see FigJT]). 
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FIG. 4: (Color Online) The entanglement of the ground state 
in the subspace L = N for bosons and AL = N for fermions. 
The interaction is repulsive Coulomb potential. 
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FIG. 5: (Color Online) Entanglement of bosons with repulsive 
Coulomb interaction, (a) Black cubic: The entanglement of 
the ground state in the subspace L = N(N — 1), namely 
fc = 1. Red circle: The entanglement of the ground state in 
the subspace L — (N — if /2 for odd N and L = N(N - 2) /2 
for even N, namely k — 2. (b) Black cubic: The entanglement 
of the ground state in the subspace L — (N~N)(N+N—3)/3, 
namely k — 3. Red circle: The entanglement of the ground 
state in the subspace L = (N - N)(N + N — 4)/4, namely 
k = 4. 



IV. BOSONS IN QUADRATIC PLUS QUARTIC 
TRAP 



So far we only consider the case in which the trap po- 
tential is a harmonic potential V(r) = ^muj 2 r 2 . This 
leads to the irrelevance of the subspace ground state 
with the interaction strength Uq, because in Eq.(f2]) 
the subspace ground state is uniquely determined by 

J2i,j,k,l=o Ui,j,k,l a l a j a k a l- Consequently, the subspace 
ground state entanglement is also irrelevant of Uq. In 
this section, we consider a quadratic plus quartic trap, 
namely 
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muj r 



1 + A (-) 
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FIG. 6: (Color Online) The eigenvalues of the single-particle 
reduced density operator of the ground state in the subspace 
AL = TV for fermions with repulsive Coulomb interaction. It 
can be seen that at TV = 18, the distribution of eigenvalues 
changes qualitatively. 
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N N 

FIG. 7: (Color Online) Entanglement of fermions with re- 
pulsive Coulomb interaction, (a) Black cubic: The entangle- 
ment of the ground state in the subspace AL — TV(TV — 1), 
namely k = 1. Red circle: The entanglement of the ground 
state in the subspace AL = (TV — l) 2 /2 for odd TV and 
AL = TV(TV - 2)/2 for even TV, namely k = 2. (b) Black 
cubic: The entanglement of the ground state in the subspace 
AL = (TV - TV) (TV + TV - 3)/3, namely k = 3. Red cir- 
cle: The entanglement of the ground state in the subspace 
AL = (TV - TV)(TV + TV - 4)/4, namely k = 4. 
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FIG. 8: (Color Online) Entanglement of bosons with con- 
tact interaction, both attractive and repulsive (A = 0.005 in 
our calculation and Uq in unit fuo). (a) The single-particle 
and two-particle entanglement of the ground state in the sub- 
space L — TV = 5. (b) The single-particle and two-particle 
entanglement of the ground state in the subspace TV = 5, 
L = N(N - 1) = 20. 

can be written as 

L L 

T-Ll = eiajai - Ltel + U ^.j.JM ! ] * *) ( 3 ) 

Z=0 i,j,k,l=0 

where 

Ui,j.k,i = JJdz 1 dz2(p*(z 1 )tp*(z 2 )U{z 1 ,Z2)(pk(z2)<Pi(zi). 

The single-particle LLL wave function tpi and its eigen- 
energy ei are solved by Numerov method. The ground 
state of Eq.(|3]) is dependent on Uq, making it possi- 
ble for us to investigate the relation between entangle- 
ment and interaction. For simplicity we suppose that 
the bosons interact with each other through contact po- 
tential (either attractive or repulsive), so that Uij t k,l — 
J dz(p* (z)ip* (z)(fk(z)(pi(z) . In order to keep the LLL ap- 
proximation valid, throughout the calculation we restrict 
the range of Uq in the interval [— 0.05Sw, 0.05hiu]. 

We consider both the single-particle and two-particle 
entanglement of the ground state in the subspace L = TV 
and L — N(N — 1). From Figs|S]and|ni we can see that 
both Si and S2 grow with the increase of |Z7o|, however, 
the curves are not symmetric with respect to the point 
Uq = 0, where no entanglement exists. In FigJSl we note 
that when the interaction is negative, the single-particle 
entanglement of the ground state in the subspace L = TV 
increases with TV, meaning that attractive bosons cannot 
condense. However, when the interaction is positive, the 
single-particle entanglement of the ground state in the 
subspace L = TV decreases with TV, leading to a conden- 
sate phase in the thermodynamic limit. 



where an = \ — is the oscillator characteristic length. 
The rotating Bose and Fermi gas confined in such an 
anharmonic potential is investigated in Refs.[26l [27|. In 
the weakly- interacting limit, the subspace Hamiltonian 



V. SUMMARY 

In this paper, we consider the particle entanglement 
in rotating Bose and Fermi gases. When the parti- 
cle number TV is fixed, for both bosons and fermions, 
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FIG. 9: (Color Online) Entanglement of bosons with contact 
interaction, both attractive and repulsive (A = 0.005 in our 
calculation and Uo in unit frw). (a) The single-particle en- 
tanglement of the ground state in the subspace L — N for 
N — 5, 10, 15. One can see that when the interaction is nega- 
tive, Si increases with N while for a positive interaction, Si 
decreases with N. (b) The single-particle entanglement of the 
ground state in the subspace L — N(N — 1) for N — 3,4,5. 
The inset of (a) is the single-particle entanglement in the sub- 
space L — N when the interaction is negative with A = 0. 

through investigating the single-particle entanglement of 
the ground state in various subspaces with fixed angular 
momentum L, we find that the phenomena of particle 



localization and vortex localization can be indicated by 
the single-particle entanglement. Moreover, we study the 
single-particle entanglement of the ground state in some 
special subspaces. For the subspace L = N, through the 
comparison between the results of bosons and fermions, 
we find that the edge reconstruction of rotating fermions 
can be indicated by the single-particle entanglement. For 
the subspace L = (N - N)(N + N - k)/k, we relate the 
subspace ground state with quantum Hall state from the 
viewpoint of entanglement and construct the same trial 
wave function as that in Rcf.|8]. For bosons, different 
properties of condensate phase and vortex liquid phase is 
reflected by entanglement. At last, we study the relation 
between entanglement and interaction by introducing a 
quartic trap. 
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